EUROPEAN LABORATORY FOR PARTICLE PHYSICS 



CERN-PPE/97-094 
24th July 1997 



Spin alignment of leading K*(892)° 
mesons in hadronic Z° decays 



The OPAL Collaboration 



Abstract 

Helicity density matrix elements for inclusive K*(892)° mesons from hadronic Z° decays have 
been measured over the full range of K*° momentum using data taken with the OPAL exper- 
iment at LEP. A preference for occupation of the helicity zero state is observed at all scaled 
momentum x p values above 0.3, with the matrix element poo rising to 0.66 ± 0.11 for x p > 0.7. 
The values of the real part of the off-diagonal element pi_i are negative at large x p , with a 
weighted average value of —0.09 ± 0.03 for x p > 0.3, in agreement with new theoretical predic- 
tions based on Standard Model parameters and coherent fragmentation of the qq system from 
the Z° decay. All other helicity density matrix elements measured are consistent with zero over 
the entire x p range. The K*° fragmentation function has also been measured and the total rate 
determined to be 0.74 ± 0.02 ± 0.02 K*(892)° mesons per hadronic Z° decay. 
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1 Introduction 



A large number of studies of inclusive vector-meson production in multihadronic Z° decay have 
so far concentrated on the measurement of fragmentation functions and total inclusive rates 
(see for example [Q and for reviews and data compilations) . Relatively little has been done 
to investigate the role of meson spin in the production dynamics. The primary quarks from 
Z° decay at LEP are highly polarized, and this may play a role in determining the helicities of 
leading vector mesons. Measurements of B* mesons have shown that these are produced 
equally in all three helicity states, as would be expected in a simple statistical picture |7|,|| 
where sea quark helicities are chosen at random and the ratio of vector to pseudoscalar meson 
production is 3:1. However, in production of D* and 0(1020) mesons at large momentum 
fraction x v (= p m cson/Pbeam), where the mesons have a high probability to contain a primary 
quark, OPAL have reported || deviations from equal populations of the three helicity states. 
The helicity density matrix element poo was found to be greater than 1/3, corresponding to the 
helicity-zero state being favoured. At the same time, the data showed evidence for negative 
values of the matrix element Repi_i for both D* and 0(1020), in agreement with qualitative 
predictions of a model based on coherent (non-independent) fragmentation of the primary quark 
and antiquark ||. Recently, DELPHI [1(J have reported measurements for p(770)°, K*(892)° 
and 0(1020), confirming the OPAL observations of vector-meson spin alignment at large x p , 
but finding no evidence for non-zero values of the off-diagonal element Repi_i. In the present 
paper, the OPAL studies are extended to the K*(892)° meson (and its antiparticle) in its K + 7r~ 
(K~7T + ) decay mode, and the helicity density matrix elements are measured with high statistics 
over the entire x p range. The findings are compared with a new theory [T]J based on a Standard 
Model description of the spin structure of the reaction e + e~ — > Z° — > qq, followed by coherent 
fragmentation of the qq pair. 



2 Formalism for the helicity density matrix 



The formalism of the helicity density matrix and its measurement using the vector-meson decay 
angular distribution have been described fully in 0. The analysis is done in the vector-meson 
helicity rest frame, where the z-axis is the direction of the K*° in the overall centre-of-mass 
system (the same as the laboratory at LEP), the y-aads is the vector product of this z-axis with 
the incident e~ beam direction, and the x-axis is such as to form a right-handed coordinate 
system. The angles #h and 0h are the usual polar and azimuthal angles of the K from the 
K*° decay, measured in this frame. The helicity density matrix elements, pw, are determined 
by fitting the observed angular distributions W using: 



M/(cos# H ) 



3 
4 



Poo) + (3poo - 1) cos 2 6> H 



W(\a\) = (2/tt) [1 + 2Repi_iCos2|a| 
W(\/3\) = (2/tt) [l + 2Imp 1 _ 1 cos2|/3| 



(1) 

(2) 
(3) 



The angles a and (3 are introduced to exploit the symmetry properties of the distributions and, 
with defined over the range — 7r to it, are given by: a = |0h| — tt/2; (3 = |0 H + 7r/4| — 7r/2 for 
0h < 37r/4; and f3 — |0h — 37r/4| — n/2 for 0h > 37r/4. Two measured asymmetries give other 
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combinations of matrix elements: 

7r AT(sin20 H cos0 H > 0) - iV(sin20 H cos0 H < 0) 
e(pio-po-i) - iV(sin 26 R cos H > 0) + N(sin 2# H cos H < 0) 



(4) 



7i Af(sin 20 H sin </>h > 0) - iV(sin 20 H sin H < 0) 
m(pio-po-i) = ^= ^ s . n 2 ^ gin H > o) + jy^in 2 Q n sin H < 0) ( ) 

where N is the number of events in the given angular range. 

The element poo of the helicity density matrix gives the relative intensity of mesons in the 
helicity state, while the off-diagonal element pi_i is a measure of coherence between the 
helicity +1 and helicity —1 states. 



3 The OPAL detector and data samples 



The OPAL detector is described in [12]. For the present analysis, the most important com- 



ponents were the central tracking chambers which consist of two layers of silicon microvertex 



detectors |13|], a high-precision vertex drift chamber, a large- volume jet chamber, and a set of 
drift chambers (the ^-chambers) which measure the coordinates of tracks along the direction 
of the beam. The OPAL coordinate system is defined with the z-axis following the electron 
beam direction; the polar angle 9 is defined relative to this axis, and r and <fi are the usual 
cylindrical polar coordinates. The central chambers lie within a homogeneous axial magnetic 
field of 0.435 T. Charged particle tracking is possible over the range | cos#| < 0.98 for the full 
range of azimuthal angles. For Kn systems around the K*° mass, the mass resolution of the 
detector varies, with momentum, from about 5 to 20 MeV. The OPAL jet chamber is capable 
of measuring specific energy loss, dE/dx, with a resolution, a (dE/dx) /(dE/dx), of 3.5% for 
well-reconstructed, high-momentum tracks in multihadronic events ||14|| . 

The present analysis used the full OPAL sample of 4.3 million multihadronic Z° decays 
recorded at LEP 1 between 1990 and 1995. To correct for losses due to the acceptance and 
efficiency of the experiment and the selection procedures, and also to provide signal and back- 
ground shapes for fits to the data mass spectra, 6 million Monte Carlo events were used, which 



had been generated using JETSET 7.4 [fLq] tuned to the OPAL data [It], and processed through 



a full simulation of the experiment [17]] and the data reconstruction and analysis. 

A detailed description of the selection of hadronic Z° decay events in OPAL is given in 
For the present analysis, tracks in the selected events were required to have: a minimum 
momentum transverse to the beam direction of 150 MeV/c; a maximum momentum of 1.07 x 
Ehe&m, based on the the momentum resolution of the detector; a distance of closest approach 
to the interaction point less than 5 cm in the plane orthogonal to the beam direction, and the 
corresponding distance along the beam direction less than 40 cm; a first measured point within 
a radius of 75 cm from the vertex; and at least 20 hits available for measurement of specific 
energy loss, dE/dx. 

Kaons and pions were identified using the dE/dx measurements. For each track, a \ 2 
probability (weight) was formed for each of the stable particle hypotheses (e, p, n, K and p). A 
track was identified as a pion or a kaon if the appropriate weight was above 5% and was larger 
than the weight for each of the other stable hadron hypotheses. Between momenta of 0.8 and 
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2.0 GeV/c, the ir, K and p bands overlap in dE/dx, leading to considerable ambiguity among 
hypotheses. Therefore no tracks were identified as kaons in this momentum range although, 
since most tracks are pions, pion identification was still allowed. With these dE/dx selections, 
pions were identified with a typical efficiency of 50% and a sample purity of 90% over most 
of the momentum range. The kaon efficiency was around 45%, while the sample purity was 
typically 30% at low momentum, rising to 50% at high momentum. 

4 Measurement of the K*(892)° helicity density matrix 
and fragmentation function 

4.1 Inclusive 'K ± 7T^ Z mass spectra 

For both the real and simulated data samples, inclusive mass spectra were formed for K ± 7r =F 
and K ± 7i ± systems in ranges of scaled momentum, x p , of the 2-particle systems. Within each of 
12 x p ranges, whose limits are given in table |l|, the data were further divided, based on the K ± 
direction in the Kir helicity frame, into each of six bins of cos#h, four bins of |a|, four bins of \(3\ 
and, for asymmetry measurements, the positive and negative regions of each of sin 2#h cos 0h 
and sin2#Hsin0H- This resulted in 12 x 18 = 216 mass spectra for each of K ± Ti Zf and K 7r ± . 
To account for combinatorial backgrounds, some of which change rapidly near the K*° peak 
region due to kinematics, background-subtracted mass spectra were formed by subtracting the 
like-charge YJ^n^ spectra from those for K ± tt t . For the Monte Carlo sample, separate spectra 
were also made for the most important states contributing to K Tr^, using information on the 
origin of each track at the generator level. 

Imperfect particle identification leads to significant contributions in the K^n^ mass spectra 
from 7T + 7r _ . The p(770)° reflection is particularly important. For kinematical reasons, its shape 
and peak position vary with x p and with the reconstructed cos #h of the misidentified pion. In 
general, the p° reflection peaks well above the mass of the K*° for negative values of cos#h but 
moves under the K*° peak as cos6*h approaches +1. Figure [I], which shows the background- 
subtracted K ± 7i :¥ mass spectra in the six cos #h bins for 0.3 < x p < 0.4, illustrates this behaviour 
of the p° reflection. The tt + 7t~ products of the uj decay also peak close to the K*° and move 
with cos^h, but they form a broader distribution. Other states, such as 0(1020) — > K + K~, are 
less important and are negligible near to the K*° peak when compared to the combinatorial 
backgrounds. 

4.2 Efficiency for K*° reconstruction 

The overall reconstruction efficiency for K*° decaying via K ± 7r =F was measured separately in 
each of the 216 bins (of x p plus an angular variable) using the Monte Carlo simulation. The 
efficiency typically varied between 10% and 25% with cos#h for all x p above 0.1, while at lower 
x p it was a stronger function of cos#h, varying from 5% to 35%. The efficiency depended much 
less on the azimuthal angles, a and j3. 

Due to the track selection and particle identification criteria, and particularly to the fact 
that no kaons were identified in the overlap region of dE/dx from 0.8 to 2 GeV/c, the K*° 
efficiency was zero for some regions in the x p versus cos^h plane, particularly for x p < 0.15. 
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The limits of the x p bins, listed in table [I], were chosen to ensure that each bin of x p contained 
a measurable K*° signal for every bin of cos#h- Because the decay particles have a momentum 
of about 300 MeV/c for a K*° at rest, measurements were possible down to x p = 0. 



4.3 Measurements of the helicity density matrix 

Each of the 216 background-subtracted mass spectra was fitted using a minimum x 2 fit, for 
100 bins of mass over the range 0.675 — 1.175 GeV, to a sum of contributions (as histograms) 
from K*(892)° signal, p(770)° reflection, u(782) reflection and background. For the K*° signal, 
the shape was taken from the detector-level Monte Carlo simulation, but because the K*° is 
generated in JETSET as a simple Breit-Wigner resonance, the output of the simulation was 
reweighted as a function of mass to reproduce a P-wave Breit-Wigner with mass-dependent 
width, which is expected to represent the K*° line shape in the real data. This procedure 
automatically takes account of variations of mass resolution with momentum and decay angles. 
The shapes of the u reflections were also taken from the simulation. To obtain the shapes of 
the residual combinatorial backgrounds for the fits, the Monte Carlo spectra were treated in 
the same way as the data, and then the contributions from K*° and the p° and u reflections 
were removed. 

The shape of the p(770)° resonance in Z° decay has long been known to be severely distorted 
from a Breit-Wigner shape Hl9[ , particularly at low x p where the peak mass is shifted by some 
40 MeV. The source of the distortion is still uncertain, but is likely to be due to residual Bose- 
Einstein correlations The standard implementation of Bose-Einstein correlations in the 



JETSET program has been found to improve significantly the agreement with data ]T9|,|2"T 



The ALEPH Collaboration p2 have measured inclusive p° production by allowing, in their 



7r + 7r~ mass spectrum fits, variation of the JETSET parameters which control the strength of 
the correlations. As the standard representation of the p° shape for the present analysis, the 
fitted ALEPH values have been used to generate p° line shapes using JETSET. Since the Bose- 
Einstein effects act primarily on low-momentum particles, shapes were generated as functions 
of both x p and cos #h since for a given value of x p the momentum of a decay pion is related to 
cos#h- As a result the line shape of the p°, and consequently its reflection in K ± 7r =F , varies with 
both x p and cos^h- To obtain the shapes of the reflections in K 7T^, the decay tracks from the 
generated p° mesons were passed through a Monte Carlo model with a simplified simulation 
of the detector and the analysis selections. The performance of this simulation was verified by 
comparing its results, for a simple Breit-Wigner shape as input, with the output of the full 
detector-level JETSET simulation. Other possible p° line shapes were considered for studies of 
systematic errors, and will be discussed in section |5.1| . 

Figure [I] shows, as an example, the results of the mass spectrum fits for the six bins of 
cos#h in the x p bin 0.3 < x p < 0.4. The average value of \ 2 f° r the fits, all of which were 
acceptable, to the 72 independent mass spectra in bins of cos#h and x p was 99.2 for 96 degrees 
of freedom. The K*° intensities extracted from the fits were fully corrected for efficiency and 
branching ratio to give, in each of the 12 bins of x p , the three inclusive double differential cross 
sections for K*° production and decay, (1 / a^)d?a / dx p dy , for y = cos^h, \oz\ and \j3\, where <7h 
is the total multihadronic Z° cross section. 

The helicity density matrix elements in each x p bin were obtained by fitting these cross sec- 
tions to the appropriate form (equations (1), (2) and (3)) as functions of the angular variables. 
Figures ^| and |3] show some of the measured differential cross sections in cos #h arid | a | along 
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with the results of the fits. The two asymmetries, in sin2#HCOs0H and sin26 l Hsin0H ; were also 
measured (equations (f|) and (H)). 

As outlined in section [4.2| the K*° efficiencies were measured for bins of cos^h and 0h> each 



integrated over the range of the other, using Monte Carlo events which were generated with 
isotropic decay distributions. Since the data clearly show non-isotropic decay distributions at 
large x p , any correlations in efficiency between cos#h and 0h could therefore bias the results 
of the measurements. The fits were therefore iterated, at each stage using the fitted values 
of the density matrix elements to generate angular distributions with which to recalculate the 
efficiencies in cos #h and 0h for the next iteration. The procedure produced only small changes, 
well within the statistical errors, to the measurements. The final results quoted are those 
obtained for two iterations, after which no further changes were induced within the quoted 
number of significant figures. 

The measured values of Im pi_i, Re (p 10 — Po-i) and Im (p 10 — Po-i) were found, as expected, 
to be consistent with zero over the entire x p range, with weighted average values for x p > 0.3 
of -0.03 ± 0.03, -0.06 ± 0.05 and -0.02 ± 0.04 respectively. For x p < 0.3 the values of p 00 
and Re pi_i were also consistent with an isotropic decay distribution. However in the large 
momentum region, x p > 0.3, the diagonal element poo was found to be larger than 1/3 while 
Repi_i was negative. Table |I| gives the measured values of p o and Repi_i in the bins of x p . 
Systematic errors, given in the table, are discussed below in section |5]. The data are also shown 
in figure f| where the clear deviations from poo= 1/3 and Repi_i= are seen for x p > 0.3, 
indicative of a 'leading particle effect', i.e. the probability for the meson to contain one of the 
primary quarks rises as the scaled momentum increases; this is discussed further in section 6. 



4.4 Fragmentation function and total rate 

For each of the 12 x p bins, the fragmentation function, (1/ 'cruder /dx p , was measured by inte- 
grating the fitted curves obtained from the cos #h spectra. The results are given in table 
The total inclusive rate for K*° production was obtained by integrating the fragmentation func- 
tion. Since the differential cross section has been measured over the entire x p range, no explicit 
systematic uncertainty arises from extrapolation using a Monte Carlo model. However there 
is an implicit assumption, borne out by the measurements, that the angular distribution in 
the Monte Carlo model is the same as that for the data for 0.03 < x p < 0.1 where there is a 
significant region of zero efficiency in the cos #h versus x p plane. The overall inclusive rate was 
determined to be 0.74 ± 0.02 ± 0.02 K*(892)° mesons per multihadronic Z° decay, in agreement 
with previous measurements |22^p4| (but with smaller statistical and systematic errors). These 



results update the previous OPAL measurements of the K*(892)° fragmentation function 
and total inclusive rate. Since the present measurements have been obtained by integrating the 
cross section, (l/crh)d 2 <j/dx p dcos^H, they are in principle more accurate at large x p than all 
previous measurements; in previous analyses the data were not binned in cos#h, and efficiencies 
were calculated assuming isotropic decay distributions. 



5 Systematic effects 
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5.1 The p(770)° line shape 



An important source of systematic error on the K*° measurements clearly arises from uncer- 
tainty in the shape of the p(770)° resonance and its reflection in the YJ^tt^ mass spectrum. 
To measure the systematic error from this source, the fits were repeated using two extreme 
possibilities for the shapes: those obtained from the simple Breit-Wigner resonance in the 
full detector-level simulation; and shapes obtained by multiplying a relativistic P-wave Breit- 
Wigner with a skewing factor which varied with x p , as used in pO[p^| . The former shape has no 
skewing of the line shape and is known to give poor fits to the ir + n~ mass spectra, particularly 
at low x p . The latter includes a shift of the peak position and a skewing of the line shape to 
low masses, but produces too sharp a cut-off at high mass. Because the p° reflection is more 
important at cos#h > 0, the entire analysis was also repeated using only the negative region 
of cos^h- For each of the measurements of inclusive cross sections and helicity density matrix 
elements the systematic error from the p° lineshape uncertainty was taken to be l/\/l2 times 
twice the maximum deviation obtained from the standard fit value. 

5.2 Detector and dE/dx simulation 

The analysis was repeated, dividing the data into two approximately equal samples covering 
two ranges of K*° production angle: | cos#| < 0.5 and | cos#| > 0.5. Here 9 is the angle between 
the K ir^ 1 system and the incident electron beam. The former range is entirely in the barrel 
region of OPAL, while the latter covers part of the barrel plus the whole of the end cap regions. 
In the barrel region, tracks may have the maximum number of jet chamber hits; at higher 
values of | cos#| the number of possible hits falls and the systematics of dE/dx measurements 
change due to the low angle of the tracks with respect to the wires of the jet chamber. As 
well as serving as a check for any systematic differences arising from the detector simulation 
of these two regions, this analysis provided results for comparison with theoretical predictions 
(see section |]). The measurements of the cross sections and p o values were in good agreement 
for the two regions, within the statistical errors, and so no systematic errors were assigned from 
this study. The values of Repi_i may be expected to vary with cos^, as will be discussed in 
section || 

To investigate further the effects of the dE / dx requirements, the minimum number of hits 
required on a track was doubled to 40, and the analysis repeated, with systematic errors taken 
to be the maximum deviation of any measurement from its standard fit value. 

Further systematic checks on the simulation of the energy loss were made by varying the 
assumed mean values of the theoretical dE/dx distributions for a given track hypothesis, and 
the assumed resolution on the energy loss measurements. Studies of well-identified pions from 
Kg decays, protons from A decays and kaons from D° and r decays were used to place limits 
on the maximum possible deviations of these quantities, and the analysis was repeated, with 
the dE/dx weights of the tracks being recalculated each time. Systematic errors were again 
assigned as 1/ \^12 times twice the maximum measured deviations from the standard fit values. 

5.3 Influence of charm decays 

Decays of charmed particles of the type P — > K*° + P', where P and P' denote any pseudoscalar 
mesons (for example D° — >K*°7r°), must produce K*° mesons which are in the helicity-zero state 
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when viewed from the P rest frame. The total contribution of all such known decays [EE] to 
the overall K*° spin alignment in the helicity frame has been studied using Monte Carlo and 
was found to be negligible. Their removal would have no measurable effect below x p = 0.2 and 
would reduce the measured values of poo by a maximum of 0.02 at x p around 0.5; their effect 
was therefore ignored. 



5.4 The K* u line shape 

To take account of possible long tails in the line shape of the K*° resonance beyond the upper 
limit of the fits to the mass spectra, the relativistic P-wave Breit-Wigner used for the fits was 
integrated out to ten full widths above the nominal peak position, resulting in an increase of 
5% in the total intensity. Since the shape of such a resonance is uncertain so far from the pole 
position, one half of this extra contribution was added to each measurement of the differential 
cross section (the results in table 1 already include these corrections), and an additional system- 
atic error of 5%/\/l2 was assigned to each measurement and to the total integrated rate. This 
systematic error has no influence on the measurements of the spin density matrix elements. 

The final systematic errors, as given in table [3], were obtained by combining in quadrature 
the contributions from all of the above sources. For all of the measurements (of cross sections 
and spin density matrix elements) the systematic errors from detector and dE/dx effects were 
found to be approximately equal to those from the p° line shape uncertainties. 



6 Comparison of results with theory 



In a recent paper [II] , a relation is derived between the values of poo and Repi_i for vector 
mesons containing a primary quark, using Standard Model parameters and some plausible 
assumptions about hadronization. According to this theory, negative values of the off-diagonal 
element Re pi_i are generated for leading vector mesons by coherence in the production and 
fragmentation of the primary quark and antiquark from the Z° decay. For the primary K*° 
mesons, it is predicted that 

sin 2 9 , , 

Repi _ lK _ . 17(1 _ Poo) __ (6) 

where 9 is the production angle of the leading K*° relative to the electron beam. In the present 
analysis, there is no explicit cut on | cos#|, although the efficiency falls rapidly for values above 
0.9. Based on the distribution of high-x p K*° production angles in the detector-level Monte 
Carlo events, the ratio Repi_i/(1 — poo) is predicted by equation (|) to be approximately 
—0.10 for mesons containing a primary quark. Table |2| shows the measured values at large x p , 
with statistical and systematic errors combined. The results in each independent x p bin are 
compatible with the theory, although there is a suggestion that the magnitudes of the measured 
values are larger than predicted. The weighted average value is —0.19 ± 0.05 for x p > 0.3. It 
should be noted however that comparison with the theory is somewhat obscured by lack of 
knowledge of the proportion of K*° mesons which contain a primary quark at a given value of 
x p . Studies using the JETSET model suggest that this fraction is around 20% at x v = 0.3, 
rising approximately linearly to 100% at x p = 1. 
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As discussed in section 5^ the density matrix elements were measured separately for two 
ranges of K*° production angle: | cos#| < 0.5 and | cos#| > 0.5. The results are given in table || 
The dependence of Repi_i on 9 predicted by equation (|6|) would imply that the ratio of these 
two measurements should be 2.3, assuming no variation of poo with cos 9. The measured ratio 
of the weighted mean values is 1.5 ± 0.7 over the range x p > 0.3, consistent within relatively 
large errors both with the theory and with no variation of Rep^! with cos#. 

While the theory of [[Oj] derives a relationship between the values of p o an d Repi_i, it does 
not offer any explanation for deviations of poo from 1/3. A number of models relevant to this 
topic are discussed in || and |l(J. In the simple statistical model based on spin counting @, 



the maximum value of poo is 0.5, and this can only occur when there is zero probability for 
production of pseudoscalar mesons containing the primary quark. While this model fits with 
the measurements of B* mesons 0-§f , it is now clearly ruled out by the LEP data for the lighter 
vector mesons. The QCD-inspired model of f26f , which predicts poo= for leading vector mesons 



is also ruled out by the measurements; however does point out that this result would only be 
expected in the limit that quark and meson masses and transverse momenta can be neglected. 
The model of ]7J, in which vector-meson production is considered as arising from the helicity- 
conserving process q— >qV, predicts poo= 1 for the leading vector mesons. Although its result 
is in accord with the observations, this model is based on a picture of parton hadronization 
which has been superseded by the more firmly based, detailed and successful string and cluster 
models; these latter models do not however address the dynamical role of meson spin in the 
hadronization. 



7 Conclusions 

From the present analysis, together with OPAL measurements of inclusive <ft and D* mesons, 
and DELPHI [|l(|] measurements of p°, K*° and 0, it is now clearly established that leading 
vector mesons, other than the B*, are spin-aligned with a preference for the helicity-zero state. 

The results of the present K*° analysis are consistent with the idea that coherence in the 
qq production and fragmentation plays a role in generating non-zero values of the off-diagonal 
element Repi_i, in agreement with the theory of [|ll|]. However the measured value of Re pi_i 
for the K*° is slightly larger than the theory predicts, while for the D* it is smaller || . This may 
indicate a flavour dependence of the mechanism producing the spin alignment, or the existence 
of more than one mechanism. While OPAL observe non-zero (and negative) values of Re pi_i 
in three separate analyses, of K*°, <fr and D* production, the DELPHI results on p°, K*° and <j) 
are all consistent with zero, within relatively small errors. Further measurements are required 
to clear up this apparent inconsistency. 



is 



The LUND string model > as implemented in the JETSET Monte Carlo program |L5 
highly successful in describing many features of hadronic Z° decay events GHJ . However, in the 



simplest string model, no spin alignment of vector mesons is expected ||27|| . Similarly, the QCD 



cluster model of HERWIG [25] has no mechanism to produce spin-aligned vector mesons; and 
neither model will generate non-zero values for the off-diagonal elements of the helicity density 
matrix. The vector meson spin alignment results therefore provide a challenge to the standard 
Monte Carlo models of the physics of inclusive hadron production in Z° decay at LEP. Can 
they be extended to include a physical mechanism which explains the measured vector-meson 
helicity density matrices? 
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Table 1: Inclusive K*° cross sections and helicity density matrix elements measured in the 
helicity-beam frame over the full range of K*° production angles. The first errors are statistical 
and the second systematic. 



x p range 


Repi_i/(l-poo) 


0.3-0.4 


-0.24 ±0.08 


0.4-0.5 


-0.17±0.13 


0.5-0.7 


-0.11 ±0.09 


0.7-1.0 


-0.24 ±0.22 


0.3-1.0 


-0.19 ±0.05 



Table 2: Ratios of K*° helicity density matrix elements measured over the full range of produc- 
tion angles at large x p . The last row gives the weighted mean value over the range indicated. 



x p range 


cos# < 0.5 


cos# > 0.5 


0.3-0.4 
0.4-0.5 
0.5-0.7 
0.7-1.0 


-0.12 ±0.04 
-0.14 ±0.06 
-0.08 ±0.05 
-0.17±0.10 


-0.13 ±0.04 
-0.05 ±0.06 
-0.02 ±0.06 
-0.01 ±0.11 


0.3-1.0 


-0.12 ±0.03 


-0.08 ±0.03 



Table 3: Values of Re pi_i measured over different ranges of K production angles at large x p . 
The last row gives the weighted mean values over the range indicated. Statistical errors only 
are given; the systematic errors largely cancel when the ratios are taken. 
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Figure 1: Inclusive background-subtracted K ± 7r :f mass spectra for the x p range 0.3 < x p < 0.4 
for the 6 bins of cos#h- The points with error bars show the OPAL data. The full histograms 
show the fits described in section 4.3, with the dashed histograms giving the K*° component 
and the shaded areas showing the contributions of the p° reflection. 
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OPAL 




Figure 2: Differential cross sections, (l/(Xh)d 2 cr/da;pdcos#H, for 6 ranges of x p . The points are 
the OPAL data and the curves show the fits used to obtain the helicity density matrix elements 
poo- The errors are statistical only. 
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Figure 3: Differential cross sections, (l/ah)d 2 a/da;pd|a:|, for 6 ranges of x p . The points are the 
OPAL data and the curves show the fits used to obtain the helicity density matrix elements 
Repi_i. The errors are statistical only. 
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Figure 4: Measured K*° helicity density matrix elements, poo an d Repi_i, as functions of x p . 
The error bars are statistical and systematic combined in quadrature, with the tick marks 
giving the statistical errors only. 
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